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Abstract 

In this paper we investigate dynamics of church bells, characterize their most common working regimes and investi¬ 
gate how to obtain them. To simulate the behavior of the yoke-bell-clapper system we use experimentally validated 
hybrid dynamical model developed basing on the detailed measurements of the biggest bell in the Cathedral Basil¬ 
ica of St Stanislaus Kostka, Lodz, Poland. We introduce two parameters that describes the yoke design and the 
propulsion mechanism and analyze their influence on the systems’ dynamics. We develop two-parameter diagrams 
that allow to asses conditions that ensures proper and smooth operation of the bell. Similar charts can be calculated 
for any existing or non-existing bell and used when designing its mounting and propulsion. Moreover, we propose 
simple and universal launching procedure that allows to decrease the time that is needed to reach given attractor. 
Presented results are robust and indicate methods to increase the chance that the instrument will operate properly 
and reliably regardless of changes in working conditions. 

Keywords: Bells, dynamics, impacting system, hybrid system, ringing schemes 


1. Introduction 


Bells are one of the oldest musical instruments which still play an important cultural role. They were invented 
in China but today their sound announces major events all around the world. Depending on the region bells 
are mounted in a number of different ways basing on local customs and tradition. In Europe we can encounter 
three characteristic mounting layouts: Central European, English and Spanish. In Central Europe bells usually tilt 
on their axis with maximum amplitude of oscillations below 90 degrees. In the English system the amplitude of 
oscillations is greater and bells perform nearly a complete rotations in both directions. Conversely, in the Spanish 
system bells rotate continuously in the same direction. All these mounting layouts were developed throughout 
centuries basing on experience and intuition of bell-founders and craftsmen. It is common that the bells are cast 
using casting mould passed down for ages from father to son and so forth. Although the design of a bell, its yoke, 
clapper and a belfry has been being improved for ages, proper modeling of their dynamics is still a challenging task. 

The dynamics a yoke-bell-clapper system is extremely complex and difficult to analyze due to nonlinear char¬ 
acteristic, repetitive impacts and complicated excitation. Moreover, depending on the mounting layout different 
dynamical states can be observed and each type of yoke has its own specific properties. However, already in 19th 
century Wilhelm Veltmann tried to describe mathematically the behavior of the famous Emperor’s bell in the 
Cologne Cathedral Em. He used simplified equations of motion and explained why the clapper always remained 
on the middle axis of the bell instead of striking its shell. His model was developed basing on the equations of a 
double physical pendulum. Heyman and Threlfrall use similar model to estimate inertia forces induced by a 
swinging bell. The knowledge of loads induced by ringing bells is crucial during the design and restoration processes 
of belfries. That is why there is a number of works considering dynamic interactions between bell towers and bells 
mounted in different manners. Muller (l2|, Steiner [l3and Schutz 0 focus on Central European mounting system 
which is also considered in the German DIN standard [l|. We can find similar studies concerning the Spanish system 
Hi and the English system 0,0. Ivorra et. al. present the comparison between the three mounting layouts and 
prove that in the Spanish system forces transmitted to the supporting structure are significantly lower than in the 
other two. Results presented in 0,0 show that in many cases we can improve bells’ working conditions slightly 
modifying the yoke or its support. 
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Apart from the studies concerning interactions between bells and their supports there is a number of publications 
focusing on the dynamics of the instrument itself or the clapper to the bell collisions. Klemenc et. al. contributed 
with a series of papers [gj, [g] devoted to the analysis of the clapper-to-bell impacts. Authors investigate the 
consequences of the repetitive hits and compare experimental data with numerical results obtained from the finite- 
element model. Presented results prove that full-scale finite-element model is able to reproduce the effect of collisions 
but requires long computational times and complex, detailed models. Therefore, it would be difficult to use such 
models to analyze the dynamics of bells. Because of that, recently we observe the tendency to use hybrid dynamical 
models which are much simpler and give accurate results with less modeling and computational effort. In dH 
authors propose lumped parameter model of the bells mounted in Central European system and prove that with 
the model we are able to predict impact acceleration and bell’s period of motion. 

In our previous publication fl3 j we present an improved hybrid dynamical model of the yoke-bell-clapper system. 
All parameter values involved in the model have been determined basing on the measurements of the biggest bell 
in the Cathedral Basilica of St Stanislaus Kostka, Lodz, Poland. Proposed model is validated by comparing the 
results of numerical simulations with experimental data. The presented results show that the introduced model is 
a reliable predictive tool and can be used for further case studies 

In this paper we describe the most common working regimes of bells and investigate how the yoke design and 
the propulsion mechanism influence their dynamics. To simulate the behavior of the yoke-bell-clapper system we 
use the hybrid dynamical model that we present in detail in our previous publication [13]. We characterize the 
solutions that can be considered as the proper operation of the bell and analyze how the geometry of the yoke and 
the driving motor output affect the dynamics of the system. We introduce two influencing parameters and develop 
diagrams that allow to asses how the bell’s behavior depends on the yoke type and propulsion amplitude. Such plots 
describe how presumed working regimes can be obtained and can be beneficial during the design and/or restoration 
processes of the bells. In addition, we investigate the time that is needed to reach given attractor. Presented results 
prove that in some cases special launching procedure of the instrument should be introduced to shorten the time of 
transient motion. Finally, we propose simple and universal control of the driving mechanism that allows to decrease 
transient time significantly. 

The paper is organized as follows. In Section 2 we describe the hybrid dynamical model of the church bell 
and introduce parameters that influence the system’s dynamics. In Section 3 we characterize the 7 most common 
working regimes and in Section 4 investigate how they can be obtained. Finally, in Section 5 the conclusions are 
drawn. 


2. Model of the system 

The hybrid dynamical model of the voke-bell-clapper system that we consider in this paper has been described 
in detail in our previous publication [131 ]. To develop the model and determine its parameters values we have 
performed detailed measurements of the existing bell named “The Heart of Lodz” (the biggest bell in the Cathedral 
Basilica of St Stanislaus Kostka in Lodz). Using the same bell we have tuned and validated the model by comparing 
the results of numerical simulations with the data collected during a series of experiments. Presented results show 
that the model is a reliable predictive tool which can be used to simulate the behavior of a wide range of yoke-bell- 
clapper systems. In next subsections we briefly describe the derivation of the equations of motion and present the 
influencing parameters. 

2.1. Geometry of the yoke-bell-clapper system 

The model that we use is build up based on the analogy between freely swinging bell and the motion of the 
equivalent double physical pendulum. The first pendulum has fixed axis of rotation and models the yoke together 
with the bell that is mounted on it. The second pendulum is attached to the first one and imitates the clapper. The 
photo of the bell that has been measured to obtain the parameters values is presented in Fig l2.ll (a). In Figs l2.ll 
(b,c) we show schematic model of the bell indicating the position of the rotation axes of the bell Oi, the clapper 02 
and presenting parameters involved in the model. For simplicity, henceforth, we use term “bell” with respect to the 
combination of the bell and it’s yoke which we treat as one solid element. 
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Figure 2.1: “The Heart of Lodz” the biggest bell in the Cathedral Basilica of St Stanislaus Ivostka (a) and its schematic model (b,c,d) 
along with physical and geometrical quantities involved in the mathematical model of the system. 


The model involves eight physical parameters. Parameter L describes the distance between the rotation axis of 
the bell and its center of gravity (point C b ), l is the distance between the rotation axis of the clapper and its center 
of gravity (point C c ). The distance between the bell’s and the clapper’s axes of rotation is given by parameter l c . 
The mass of the bell is described by parameter M, while parameter Bb characterizes the bell’s moment of inertia 
referred to its axis of rotation. Similarly, parameter m describes the mass of the clapper and B c stands for the 
clapper’s moment of inertia referred to its axis of rotation. 

Considered model has two degrees of freedom. In Fig. 12.11 (d) we present two generalized coordinates that we 
use to describe the state of the system: the angle between the bell’s axis and the downward vertical is given by 
ipi and the angle between the clapper’s axis and downward vertical by ip 2 - Parameter a (see 12.11 (d)) is used to 
describe the clapper to the bell impact condition which is as follows: 

|yi-H = a (2-1) 

Synonymously, collision between the bell and the clapper occurs when the absolute difference between the bell’s 
and the clapper’s angular displacements is equal to a. 


2.2. Equations of motion - modeling of an oscillatory motion of the system 

In this section we present the mathematical model that we use to simulate oscillatory motion of the investigated 
yoke-bell-clapper system. We use Lagrange equations of the second type and derive two coupled second order ODEs 
that describe the motion of the considered system (full derivation can be found in 0): 


(.Bb + ml 2 ) ipi + ml c lip2 cos {ip2 — ipi) — ml c lip 2 sin (ip2 — ipi) + (ML + ml c ) gsintpi 
+D b ip i - D c (ip 2 - <Pi) = M t (ipi), 


( 2 . 2 ) 


B c ip 2 + ml c lipi cos (ip 2 ~ fi) + ml c lip 2 sin (ip 2 — ipi) + mgl sin ip 2 + D c (ip 2 — ipi) = 0. (2-3) 

where g stands for gravity and M t (ipi ) describes the effects of the linear motor propulsion. The motor is active - 
and excites the bell - when its deflection from vertical position is smaller than 7 t/15 [rad] (12°). The generalized 
momentum generated by the motor M t (ipi) is given by the piecewise formula: 


[ T sgn(y>i) cos (7.5y>i), if \ipi\ < ^ 

Mt(ipi) = l (2.4) 

l if l^i I > te 

where T is the maximum achieved torque. Although, the above expression is not an accurate description of the 
effects generated by the linear motor, in [l3j we prove that it is able to reproduce the characteristics of the modern 
bells’ propulsions. 

There are eleven parameters involved in the mathematical model presented above. The parameters have the 
following values: M = 2633 [kg], m = 57.4 [kg], Bb = 1375 [kgm 2 ], B c = 45.15 [kgm 2 ], L = 0.236 [m], l = 0.739 [m], 
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l c = —0.1 [to] and a = 30.65° = 0.5349 [rad], D c = 4.539 [Nms], D b = 26.68[IVms], T = 229.6[Nm]. As afore¬ 
mentioned, all parameters values have been evaluated specifically for the purpose. For integration of the model 
described above we use the fourth-order Runge-Kutta method. ODEs [2721 and [2731 together with the discreet model 
of impact described in the next Subsection create a hybrid dynamical system. 


2.3. Modeling of the clapper to the bell impact 

After considering free motion conditions, in this Section we briefly refer to the discreet impact model which is 
was proposed by Meneghetti and Rossi [lH and validated experimentaly in our previous publication 0- When the 
condition no is fulfilled we stop the integration process. Then, instead of analyzing the collision course, we restart 
simulation updating the initial conditions of equations 12.21 and 12.31 by switching the bell’s and the clapper’s angular 
velocities from the values before the impact to the ones after the impact. The angular velocities after the impact 
are obtained taking into account the energy dissipation and the conservation of the system’s angular momentum 
that are expressed by the following formulas: 


1 2 1 2 

-B c {tp 2 ,Ai - <Pi,ai) = k^B c (ip2,BI ~ <Pi,bi) ■ 


(2.5) 


[B b + mil + ml d cos (ip 2 - ¥>i)] ipi'Bi + [Be + ml c l cos (ip 2 - Vi)\ T2 ,bi = 
[B b + mil + ml d cos (ip 2 - Vi)\ <Pi,ai + [B c + ml c l cos (ip 2 ~ V\)\ f>2 ,ai 


( 2 . 6 ) 


where index AI stands for “after impact”, index BI for “before impact” and parameter k is the coefficient of energy 
restitution. In our simulations we assume k = 0.05 referring to a series of experiments performed by Rupp et. al. 
(l5^ . In our previous investigation [13] we have analyzed influence of k on the response of system and we have 
proved that dynamics of the system barely changes for small alterations of parameter k (±20%). Hence, we claim 
that there is no need to further adjust the value of k for the considered bell. 


2.f. Influencing parameters 

The mathematical model of the ringing bell described in Section [2] contains 12 parameters ( B b , M, L, D b , a, 
B c , to, l, D c , l c , k, T ) but most of them are self dependent. Moreover, we have to remember that we consider a 
musical instrument, hence we cannot change some of its features as it could affect the sound it generates. In real 
applications we can easily modify the driving motor and the mounting of the bell (by changing the design of the 
yoke). Therefore, we try to describe how changes of the propulsion and the yoke’s design influence the system’s 
dynamics. As a reference, we use values of parameters characteristic for “The Heart of Lodz” and alter them to 
simulate the changes in the propulsion or mounting layout. 

In our investigation we assume the linear motor propulsion that is described by piecewise function given 

by the formula 12.41 Hence, we can modify the effects generated by the motor by changing the range of bell’s 
deflection in which the motor is active (in our case (—7r/15, 7r/15)) or by altering the maximum generated torque 
T. Practically, it is much easier to modify the maximal momentum which can be realized either by replacing the 
motor or by changing the length of the torque arm. Therefore, we take T as the parameter that describes the 
driving motor characteristic and use it as the first controlling parameter. 

To describe the modifications of the yoke we introduce new parameter l r and take it as the second influencing 
parameter. In Fig. 12.21 we explain the meaning of l r parameter. 

For the reference yoke (“The Heart of Lodz” yoke) we assume l r = 0. If we change the yoke design so that the 
bell’s center of gravity is lowered, then l r < 0 and as the value of l r we take the distance by which the bell’s center 
of gravity is shifted with respect to the reference yoke. Similarly, if the bell’s center of gravity is elevated, we assume 
l r > 0 and take its displacement as the value of l r . As a maximum considered value of l r we take 0.235 [to] because 
for l r = 0.236 [to] rotation axis of the bell goes through its center of mass. 

It is important to point out that changes of l r value affect other parameters of the model. Hence, whenever we 
change the value of l r we have to swap three parameters: L has to be replaced by L r , l c must be replaced by l cr , 
and finally B b by B br . Values of L r , l cr and B br are calculated using the following formulas: 

Lj j. — L ly~ 


t C r L c L r 

B br = (B b - ML 2 ) + ML 2 r 


(2.7) 
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(a) (b) (c) 



Figure 2.2: The meaning of parameter l r : (a) reference yoke design l r = 0 (based on “The Heart of Lodz” yoke), (b) l r < 0, (c) l r > 0. 


3. Different working regimes of the bell 

The considered system is not only hybrid and nonlinear but also piecewise due to the model of driving linear 
motor. Hence it can exhibit a plethora of different dynamical phenomena such as periodic, quasi periodic and chaotic 
attractors with different number of collisions. Therefore, a particular bell can behave quite differently depending 
on the yoke design and the driving technique. Still, only few types of periodic attractors can be considered as a 
proper working regimes and have practical applications. These regimes are often called ringing schemes and can be 
classified in groups that have common characteristics. The features that are important during the categorization 
are principally: the period between the consecutive impacts, the number of collisions during one period of motion 
and the course of the collisions. In this section we point out and describe the seven most commonly meet types 
of behavior. We use four types of plots to present considered working regimes: phase portraits of the bell and the 
clapper, trajectory projected in the section of phase space showing relation between the angular position of the bell 
and the clapper and time trace of the clapper’s velocity. We use blue lines to present attractors and trajectories; 
red lines and arrows are used to present changes of velocity that are the effects of collisions while green lines are 
plots of the impact condition 12.11 and indicate when the collisions occur. 

3.1. Falling clapper 

We say that the bell works in a “falling clapper” manner if the collisions between the bell and the clapper occur 
when they perform an anti-phase motion. This type of behavior is common for bells that are mounted in the 
European manner and “The Heart of Lodz” is an example of a bell that works in this regime. In the falling clapper 
ringing scheme the amplitude of the clapper’s motion is smaller than the bell’s.The clapper’s velocity sign changes 
when collision occurs. The amount of energy that is transferred during the impact is relatively large and sufficient 
to let the bell resound nicely. Thanks to this feature it is relatively easy to receive nice voicing of the bells working 
in falling clapper regime. But, in some cases, the energy transfer is too abrupt which can lead to crack or damage 
of the bell or the clapper. 



Figure 3.1: Presentation of the symmetrical “falling clapper” with 2 impacts per period obtained for T = 350 [:Vm] and l, = 0.05 [m]. 
Phase portraits of the bell (a), the clapper (b), section of phase space with green lines that indicate impact conditions (c) and time 
trace of the clapper’s velocity (d). Effects of collisions (changes of velocity) are marked with red lines. 
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We can distinguish a symmetric type of falling clapper with 2 collisions per one period of motion and its 
asymmetric version with 1 impact per period. Although the course of the collision is similar in both cases they 
are completely different for the listeners as the time interval between successive impacts is almost doubled in the 
asymmetric case. In Figs . 13.II and 13.21 we present two types of falling clapper ringing scheme. 


(a) (b) 



(c) (d) 



<Pi [rad\ <[s] 


Figure 3.2: Presentation of the asymmetrical “falling clapper” with 1 impact per period obtained for T = 150 [IVm] and l r = —0.03 [m]. 
Phase portraits of the bell (a), the clapper (b), section of phase space with green lines that indicate impact conditions (c) and time 
trace of the clapper’s velocity (d). Effects of collisions (changes of velocity) are marked with red lines. 


3.2. Flying clapper 

If the collisions between the bell and the clapper occur when they perform in-phase motion we say that the bell 
works in a “flying clapper” manner. In this regime the amplitude of the clapper’s motion is larger than the bell’s 
and the clapper’s velocity sign remains the same after the collisions. The collisions have more gentle characteristic 
and the amount of energy that is transferred during the impact is much smaller than in “falling clapper” ringing 
scheme. Therefore, sometimes it may be difficult to achieve nice resounding of the bell, but the risk of the crack or 
damage is significantly decreased. In Figs. I3.3l and [3.4l we present two types of “flying clapper” behavior: symmetric 
attractor with 2 impacts per period and asymmetric one with only 1 impact per period respectively. 



(p x [ rad] (p 2 [rad] (p x [rad] t [5] 


Figure 3.3: Presentation of the symmetrical “flying clapper” with 2 impacts per period obtained for T = 450 [Nm] and l r = —0.91 [m]. 
Phase portraits of the bell (a), the clapper (b), section of phase space with green lines that indicate impact conditions (c) and time 
trace of the clapper’s velocity (d). Effects of collisions (changes of velocity) are marked with red lines. 


3.3. Double kiss 

The “double kiss” is the name of working regime in which we observe 4 impacts per one period of motion. 
During one period the clapper hits each side of the bell’s shell twice. The first collision on each side is in the “falling 
clapper” manner while the second impact has “flying clapper” course. This behavior is especially attractive for the 
listeners. Still, it is rarely met mainly because it is difficult to achieve. In the next Section we precisely determine 
conditions (ranges of parameters T and l r ) for which we can observe this working regime. In Fig. 13.51 we show the 
characteristics of “double kiss”. 
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Figure 3.4: Presentation of the asymmetrical “flying clapper” with 1 impact per period obtained for T = 325 [Aim] and l r = —1.21 [m]. 
Phase portraits of the bell (a), the clapper (b), section of phase space with green lines that indicate impact conditions (c) and time 
trace of the clapper’s velocity (d). Effects of collisions (changes of velocity) are marked with red lines. 



Figure 3.5: Presentation of the “double kiss” phenomenon with 4 impacts per period. The plots were obtained for T = 175 [Aim] and 
l r = 0.16 [m[. Phase portraits of the bell (a), the clapper (b), section of phase space with green lines that indicate impact conditions 
(c) and time trace of the clapper’s velocity (d). Effects of collisions (changes of velocity) are marked with red lines. 


3-4- Sticking clapper - sliding dynamics 

“Sticking clapper” is the name of working regime in which the clapper and the bell remain in contact for a 
certain amount of time. In other words the “sticking clapper” refers to the attractors that contain sliding mode. 
This working regime is typical for the bells mounted and operated in the English manner. In the considered system 
prior the sliding mode we observe a number of successive impacts (usually 3) that have a “falling clapper” course. 
The energy amount that is transferred between the bell and the clapper decreases with each subsequent collision. 
Hence, the sound effects caused by each hit are different and not all collisions may be noticed by the listener. 
Moreover, when the clapper remains in the contact with the bell’s shell it also influences the produced sound. All 
these features have to be taken into consideration if we want the bell to sound nicely. 



Figure 3.6: Presentation of the “sticking clapper” phenomemon. The plots were obtained for T = 125 [A/m] and l r = 0.2 [m[. Phase 
portraits of the bell (a), the clapper (b), section of phase space with red lines that indicate impact conditions (c) and time trace of the 
clapper’s velocity (d). Effects of collisions (changes of velocity) are marked with red lines. 


3.5. No impacting attractors and other possible solutions 

If the forcing amplitude is not sufficient or the yoke is designed improperly, we can observe stable periodic 
attractor with no collisions. In such conditions no sound is produced and the bell can not work as a musical 
instrument. Unfortunately, no impacting attractors occur in wide range of T and l r and even nowadays it is 
common that we have to redesign the yoke to make ringing possible. Apart from the ringing schemes described 
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in this Section we can indicate other periodic attractors that can be successfully employed. For example we can 
observe asymmetric flying clapper behavior with doubled period and 4 impacts per period. Such behavior can 
be easily taken as a typical falling clapper. Similarly, a quasi-periodic attractor with almost equal time intervals 
between the subsequent impacts can sound almost like a periodic ringing scheme. Moreover, as we proved in our 
previous publication 0, the course of successive collisions can differ a bit. Hence, in practical applications we 
do not demand strictly periodic behavior. Instead, we want to receive similar collisions of a presumed type in a 
possibly equal time intervals. Still, when designing the yoke and the propulsion of the bell we should always look 
for conditions that ensure periodic behavior as it increases the probability that the real system will work properly 
and reliably. 


4. Influence of the yoke design and forcing amplitude on the system’s dynamics 

4-1. One parameter diagram - transitions between different working regimes 

In Fig. 14.11 we present one parameter diagram that presents how the system behaves with varying value of l r 
parameter. We assume T = 412.5 [Nm] and analyze the dynamics of the system for l r e (—1.3, 0.23} [m\. Each time 
we start integration form zero initial conditions (tpi = 0, y> 2 = 0, ipi = 0, </j 2 = 0), hence it is not a bifurcation 
diagram (the bell and the clapper always start their motion from hanging down position with zero velocities). We 
investigate the behavior of the system after 38 minutes of preliminary motion (the analysis of influence of transient 
time on the system’s dynamics is shown in the end of this Section). To determine in which regime the system 
operates, we use two indicators. The first is the absolute value of the clapper’s angular velocity |<^ 2 | when the 
clapper passes the downward vertical <p 2 = 0. The second is the clapper’s velocity just after the clapper to the 
bell impact dpi.Ai- Moreover, we calculate the energy dissipated at instant of collision Ed- We record these values 
during the period of two minutes. Subplots (a) and (b) of Fig. 14. II correspond to the first and the second indicator 
respectively while in lower part of the Figure (subplot (d) ) we present the energy dissipation. Analyzing Fig. 14.11 
(a,b) we can asses the periodicity and the symmetry of the attractor. Moreover, we are able to evaluate the number 
of collisions per period of motion. To determine if we observe the “flying clapper” or “falling clapper” ringing scheme, 
we have to look at the phase portrait of the clapper or at least value of the clapper angular velocity just before the 
impact ip 2 ,bi- We focus on the clapper’s motion because the changes in the bell’s behavior are often imperceptible 
due to its much bigger mass and inertia. 

The two most common ringing schemes are symmetric “flying” and “falling” clapper. Analysing subplot (c) of 
Fig. 14.11 we see that in “flying clapper” working scheme collisions have gentle course and the risk of damage of the 
bell or the clapper is minimal. In this regime the amount of energy dissipated during collision is relatively small and 
do not depend on the yoke design. Conversely in “falling clapper” ringing scheme the amount of energy dissipated 
at instant of collision strongly depends on the yoke design. Collisions are more abrupt but still the bell can operate 
reliably if the l r and T parameters are properly chosen. The risk of failure is also higher for asymmetric attractors 
(both “falling” and “flying” clapper) because then the amount of energy transferred during collisions is significantly 
higher than in symmetric case (for similar parameters’ values). That is why these working regimes are rarely met. 

Vertical lines in Fig. 14.11 mark the ranges of l r where given working regimes can be observed. We see that 
transitions between different states of the system are sudden and sometimes even minor changes in the yoke geometry 
can result in a completely different behavior of the system. We usually expect the bells to work reliably for decades 
without any maintenance. Therefore, the yoke and the propulsion of the bell should be designed to ensure that the 
instrument will work properly regardless of small changes of influencing parameters. Such approach requires the 
in-depth knowledge of the systems dynamics. 

4-2. Two parameters ringing scheme diagrams 

To examine how the yoke’s design - described by parameter l r - and the amplitude of forcing - parameter T - 
influence the system’s dynamics we perform the series of numerical simulations. In our analysis we consider the 
following ranges of these parameters: l r ef- 1.3, 0.23) [m] and Te(100, 625) [Nm]. We take 154 equally spaced 
values of l r (with the step 0.01 [m]) and combine each with 42 values of T (from 100 [Nm] to 625 [JVm] with the 
step 12.5[lVm]). That gives us 6468 different sets of system’s parameters. For each one we simulate the system’s 
motion starting from zero initial conditions (ip i = 0, <^> 2 = 0, tfii =0, <p 2 = 0). 

In Figs. 14.31 and 14.21 we show the possible behavior of the systems 2 minute and 38 minutes after the excitation 
started respectively. We call these two parameters plots as the “ringing scheme diagrams” because they allow to 
determine which ringing scheme can be achieved for the given values of T and l r . We consider 7 most characteristic 
types of the bells behavior that are described in detail in the previous Section. Of course many different stable 
attractors can exist for a given set of parameters l r and T, but during our investigation we concern only the solutions 



2 impacts per period: 

1 impact per period: 

6 - symmetric falling clapper 

2 - symmetric flying clapper 

5 - asymmetric falling clapper 

1 - asymmetric flying clapper 


3 - non periodic 

4 - no impacts 

7 - double kiss 

8 - sticking clapper 


Figure 4.1: One parameter diagram showing the behaviour of the yoke-bell-clapper system with varying value of l r parameter. Supblot 
(a) presents clapper’s angular velocity when the clapper passes the downward vertical, (b) shows clapper’s velocity just after the impact 
and subplot (c) presents the energy dissipated at instant of collision. 
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that basins of attraction contain zero initial conditions = 0, <p 2 = 0, tpi = 0, </>2 = 0). As we mention before, 
we do so because in most cases the bell and the clapper start their motion from hanging down position with zero 
velocities. Another aspect that we consider is the transient time needed to reach a given attractor. It is important 
in real applications because the launching procedure of the bell should be possibly short. Hence, when designing 
the yoke and the propulsion mechanism we have to know both how presumed ringing scheme can be achieved as 
well as the time of transient behavior before we reach that attractor. If this time is too long we can either change 
the yoke design or apply special starting procedure - by the control of the driving motor. 



2 impacts per period: 

1 impact per period: 

| symmetric falling clapper 
| symmetric flying clapper 

1 1 asymmetric falling clapper 

H asymmetric flying clapper 


| no impacts 
I I double kiss 
I I sticking clapper 


Figure 4.2: Two parameter ringing schemes diagram showing the behavior of the system 38 minutes after the excitation started. 


The first plot - Fig. 14.21 - presents the state of the system 38 minutes after the propulsion starts. It shows which 
stable attractor will be achieved starting from zero initial conditions for different values of l r and T. Hence, Fig. 
14.21 indicates which ringing scheme can be obtained with given type of yoke and maximal driving torque. Analyzing 
the plot we can draw the following conclusions. The biggest areas correspond to the three most common ringing 
schemes: “falling clapper” (black) and “flying clapper” (red) with 2 impacts per one period of motion and “sticking 
clapper” (pink) that is typical for the bells mounted in the English manner. This means that these behaviors are 
relatively easy to achieve and remain even after long period of time when values of some parameters can change a 
bit (for example maximum torque generated by the linear motor). 

Generally, the design of the yoke determines how the bell will operate. For l r < —0.4 [to] we can observe “flying 
clapper” with 2 (red) or 1 (green) impacts per period but these ringing schemes can be achieved only when T is bigger 
than some threshold. The threshold value decreases with the decrease of l r that causes the growth of the inertia 
of the yoke-bell assembly with respect to its axis of rotation. If T is not sufficient we will not observe any impacts 
(blue) or the system will reach different attractor (white) - periodic or non-periodic one. For l r e (—0.4. — 0.284) [to] 
we cannot obtain any of the analyzed ringing schemes no matter how big the forcing amplitude is. This proves that 
when the yoke is designed improperly it may be impossible to force the system to ensure proper operation. For 
l r > —0.284 [to] the system can work in the following manner: “falling clapper” with 1 (yellow) or 2 (black) impacts 
per period, “double kiss” (light blue) or “sticking clapper” (pink). Each of these ringing schemes can be achieved 
despite the value of T but the yoke should be designed for the purpose. Hence, for these working regimes there is 
no need to use very powerful driving motors and the system can work more efficiently. 
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Figure 4.3: Two parameter ringing schemes diagram showing the behavior of the system 2 minutes after the excitation started. 


In real applications we usually cannot accept launching times bigger than two minutes. To determine for which 
parameters the start up procedure runs quickly, we calculate the second plot - Fig. 14.31 - that shows the state of 
the system two minutes after the propulsion started. In this case we analyze the system’s motion forl20 [s] and 
do not demand to reach the stable solution. Instead, we check if for the average observer the behavior of the 
system could be seen as proper and regular. For example we do not care if the amplitude of the bell’s motion is 
still slightly rising and just check if the period between the successive impacts is constant and the collisions have 
recurrent course. Comparing Figs. 14.21 and 14.31 we see that the difference in the size of areas is especially visible for 
l r < —0.4 [to]. The range that corresponds to the attractor without impacts (blue color) is significantly larger and 
we cannot observe any “flying clapper” ringing schemes. This means that these working regimes not only requires 
larger amplitude of forcing but also longer transient time before the system reaches presumed solution. Therefore, 
in practical applications we often apply special launching procedure to decrease the time of transient behavior. 
We also see that for l r > —0.4 [m\ the areas corresponding to different ringing schemes are almost identical to the 
ones presented in Fig. 14.21 The main difference is that for T > 230 [Am] there is a white region between the 
areas corresponding to asymmetric (yellow) and symmetric (black) “falling clapper”. Therefore, special starting 
techniques are rarely required if the yoke design ensures small inertia with respect to the bell’s rotational axis. 

Another important conclusion from the analysis of the ringing scheme diagrams is the possibility to design the 
yoke and the propulsion to enable more than one ringing scheme. For example, assuming that l, r = 0 we will observe 
“falling clapper” with 1 impact per period for Te (100, 191) [./Vm] and 2 impacts per period for Te (191, 625) [Am]. 
Hence, if the adopted driving motor enables to adjust the output torque in a sufficiently wide range, we can change 
the way the bell operates. Thanks to that, we are able to differentiate the tune of the bell for different occasions. 

4.3. Example of the launching procedure 

The results presented in the previous Subsection prove that we often have to use special launching procedure 
in order to achieve presumed ringing schemes and maintain considerably short starting time. In this Section we 
propose a simple control of the driving motor that enables to achieve all considered working regimes after only two 
minutes. We assume that launching procedure takes 120 seconds and over that time we adjust the maximum torque 
that is generated by the motor. Hence, during the first 120 seconds, instead of using formula 12.41 to describe the 
generalized momentum generated by the propulsion mechanism, we use the following function: 
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where: t is time and Di = 200 [Nm 2 ], D 2 = 24000 [Nm 2 s] are control parameters. We obtained the above 
formula basing on the following assumptions. The amplitude of forcing decreases as a square function and for 
t = 120 [s] it is equal to parameter T. Hence, there is no jump in the torque generated by the motor when the 
launching procedure ends (t = 120 [s]) and we switch back to formula [2.41 Analyzing Figs. 14.21 and 14.31 we see that 
the launching procedure is needed especially for small l r and large T. Therefore, we assume that the maximal value 
of generated torque should increase with the increase of T and the absolute value of l r . 



Figure 4.4: Maximum torque generated by the linear motor during the first 3 minutes of motion. 


In Fig. 14.41 we show how the maximum torque generated by the linear motor (for T = 500 [Nm] and l r = 
—0.44 [m]) changes during the analyzed first 3 minutes of the system’s motion. The dotted vertical line in Fig. 14.41 
marks the time when the launching procedure ends and the output of the motor stabilizes. 
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Figure 4.5: Two parameter ringing schemes diagram showing the behavior of the system with applied launching procedure (2 minutes 
after the excitation started). 


In Fig. 14.51 we present the ringing scheme diagram obtained for the system with applied launching procedure 
that is described above. Similarly to Fig. 14.31 this time we also do not require the system to reach the stable 
attractor and just check if its behavior could be perceived as proper by an average listener. Analyzing Fig. 14.51 we 
see that even using simple control of the driving motor we can achieve satisfactory results and decrease the time of 
the transient motion. Moreover, obtained results prove that formula 14.11 - universal for all T and l r - allow us to 
achieve good results in a wide range of parameters values. But, the proposed control is not efficient with respect 
to the systems with large forcing amplitude and yoke design that provides large inertia of the bell (left upper part 
of the plot.). Comparing Figs. 14.21 and 14.51 we see the difference in the borders and sizes of the areas corresponding 
to different working regimes, but still thanks to the proposed launching procedure we can reach all 7 considered 
working regimes. Hence, the proposed control allow us to obtain all considered ringing schemes but they appear in 
smaller ranges of parameters values. Therefore, in some cases there is a need to adjust the launching procedure to 
shorten transient time and reach presumed behavior of the system. 


5. Conclusions 

Hybrid dynamical model of the yoke-bell-clapper system can exhibit a plethora of different dynamical behaviors 
but only a few of them can be considered as a proper working regimes of the instrument. Therefore, the yoke of the 
bell and its propulsion mechanism should be designed to ensure the system will behave correctly and sound nicely. 
In case of bells with typical shape and mass engineers usually design yokes basing on their experience and local 
tradition. Still, it is often a challenge to create a reliable suspension for unique bells and sometimes the mounting 
of the bell or its propulsion have to be redesigned to achieve presumed ringing scheme. 

In the paper we present a method that can be used to determine the conditions under which given type of 
behavior can be achieved. We develop two parameter ringing scheme diagrams that describe how the geometry of 
the yoke and maximum output of the driving motor influence the dynamics of the system. Similar charts can be 
calculated for any bell and used to design its mounting and propulsion. Thus, we can ensure that the instrument 
will work properly and reliably regardless of small changes of parameters. 

Moreover, analyzing the ringing scheme diagrams we see that it is possible to design the yoke and the propulsion 
to enable more than one ringing scheme with relatively small changes of motor and/or yoke parameters. Hence, we 
can easily differentiate the tune of the bell for different occasions. 
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Basing on the presented analysis we are able to determine the minimum required power of the driving motor 
and estimate the time that is needed to achieve given attractor. If the transient time is too long, special launching 
procedure has to be applied. We propose the control of the driving motor that enables to reach all considered 
ringing schemes after 120 seconds of preliminary motion. It is possible to further optimize the launching procedure, 
but to get best results we need to know the values of all system’s parameters. Still, presented results are robust 
and prove that even simple and universal formula can be highly efficient in a wide range of systems’ parameters. 

Acknowledgment 

This work is funded by the National Science Center Poland based on the decision number DEC-2013/09/N/ST8/04343. 
We would especially like to thank the Parson of Cathedral Basilica of St Stanislaus Kostka Prelate Ireneusz Kulesza 
for his support and unlimited access to the bell. We have been able to measure the bell’s template thanks to the 
bell’s founder Mr Zbigniew L. Felczynski. The data on the clapper, the yoke and the motor have been obtained 
form Mr. Pawel Szydlak. 

References 

References 

[1] Glockenti^cerme DIN 4178:2005-04, 2005. 

[2] S. Ivorra and J. R. Cervera. Analysis of the dynamic actions when bells are swinging on the bell-tower of bon- 
repos i mirambell church (valencia, Spain). IProc. of the 3rd international seminar of historical constructions, 
pages 413-19, 2001. 

[3] S. Ivorra and F. J. Pallares. Dynamic investigations on a masonry bell tower. Engineering Structures, 28(5): 

660 - 667, 2006. 

[4] S. Ivorra, F. J. Pallares, and J. M. Adam. Dynamic behaviour of a modern bell tower - a case study. Engineering 
Structures, 31(5):1085 - 1092, 2009. 

[5] BD. Threlfall J. Heyrnan. Inertia forces due to bell-ringing. International Journal of Mechanical Sciences, 18: 
161-164, 1976. 

[6] A. Selby J. Wilson. Durhamm cathedral tower vibrations during bell-ringing. In Proceedings of the conference 
engineering a cathedral., pages 77 100. London: Thomas Telford, 1993. 

[7] A. Selby J. Wilson. Dynamic behaviour of masonry church bell towers. In Proceedings of the CCMS symposium, 
pages 189-199. New York: Chicago ASCE, 1997. 

[8] J. Klemenc, A. Rupp, and M. Fajdiga. A study of the dynamics of a clapper-to-bell impact with the application 
of a simplified finite-element model. Engineering with Computers, 27(3):261-272, 2011. 

[9] J. Klemenc, A. Rupp, and M. Fajdiga. Dynamics of a clapper-to-bell impact. International Journal of Impact 
Engineering, 44(0):29 - 39, 2012. 

[10] M. Lepidi, V. Gattulli, and D. Foti. Swinging-bell resonances and their cancellation identified by dynamical 
testing in a modern bell tower. Engineering Structures, 31(7):1486 - 1500, 2009. 

[11] G. Meneglietti and B. Rossi. An analytical model based on lumped parameters for the dynamic analysis of 
church bells. Engineering Structures, 32(10):3363 - 3376, 2010. 

[12] FP. Muller. Dynamische und statische gesichtspunkte beim bau von glockenturmen. Karlssruhe: Badenia 
Verlag GmbH, pages 201-212, 1986. 

[13] P. Perlikowski P. Brzeski, T. Kapitaniak. Experimental verification of a hybrid dynamical model of the church 
bell. International Journal of Impact Engineering, 2015. 

[14] KG. Scliutz. Dynamische beanspruchung von glockenturmen. Bauingenieur, 69:211-217, 1994. 


14 


[15] R. Spielmann, A. Rupp, M. Fajdiga, and B. Aztori. Kirchenglocken-kulturgut, musikinstrumente und 
hochbeanspruchte komponenten. Glocken-Lebendige Klangzeugen. Schweizerische Eidgenossenschaft, Bunde- 
samt fur Kultur BAK., pages 22-39, 2008. 

[16] J. Steiner. Neukonstruktion und sanierung von glockenturmen nach statischen und dynamischen gesichtspunk- 
ten. Karlssruhe: Badenia Verlag GmbH , pages 213-237, 1986. 

[17] W. Veltmann. Uerber die bewgugn einer glocke. Dinglers Polytechnisches Journal , 220:481-494, 1876. 

[18] W. Veltmann. Die koelner kaiserglocke. enthullungen uber die art und weise wie der koelner dom zu einer 
mii'^oerathenen glocke gekommen ist. Hauptmann, Bonn , 1880. 


15 


